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ABSTRACT 

The observational tests for the period of reionization of the universe are discussed. 
We show that this period can be observed as negative intensity patches of the CMB 
radiation with the amplitude ST/T ~ 10~ 5 and the angular sizes 0t ~10 angular 
seconds in range of the wavelength 0.1 cm< A <1 cm. The expected number density 
and frequency dependence of the amplitude permit to recognize this effect and to 
discriminate it from the noise. This method applied to the small scale variations of 
CMB temperature complements well the traditional spectral approach. 

The number density and the amplitude of observed 'negative' intensity patches 
depend upon the redshift of reionization that allow to estimate roughly this redshift. 
The ionized bubbles formed just before the period of reionization could also be seen as 
the highest peaks. The expected results are sensitive to the Jeans scale at the period 
of reionization, to the small scale shape of the primordial power spectrum and to the 
mass of dark matter particles. 
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1 INTRODUCTION 

The analysis of the cosmic microwave background radiation 
(CMB) anisotropy generated after reionization of the uni- 
verse has been began in Sunyaev & Zel'dovich (1980) (static 
and kinematic Sunyaev-Zel'dovich effect), Kaiser (1984) and 
Ostriker & Vishniac (1986). Later on various aspects of 
this problem were discussed in many publications (see, e.g., 
Tegmark, Silk, Blanchard, 1994; Hu, Scott & Silk 1994; Persi 
1995; Tegmark & Silk 1995; Aghanim et al. 1996; Hu & 
White (1996); Knox, Scoccimarro & Dodelson 1998; Gruzi- 
nov & Hu 1998; Peebles & Juszkiewicz 1998, Griffiths, Bar- 
bosa & Liddle 1999; Haiman & Knox 1999; Peterson et al. 
1999; Weller, Battye & Albrecht 1999; Rees 1999; Molnar 
& Birkinshaw 2000; Cooray, Hu, & Tegmark 2000; Benson 
et al. 2001). Various possible sources of such perturbations 
were considered - from linear and second order perturba- 
tions and up to CMB scattering within a set of high density 
clouds and clusters of galaxies. 

Firstly the expected variations of the CMB temper- 
ature have been calculated in Vishniac (1987) where has 
been noted that, for larger I, the power spectrum of the 
temperature variations reproduces well the spectrum of per- 
turbations and is defined by the reionization epoch. Recent 
publications repeat the same calculations, concentrate more 



attention on the impact of nonlinearity and inhomogene- 
ity of reionization (Peebles & Juszkiewicz 1998; Gruzinov 
& Hu 1998; Knox, Scoccimarro & Dodelson 1998; Jaffe & 
Kamionkowski 1998; Haiman & Knox 1999; Hu 2000), and 
use numerical simulations of the reionization (Springel et 
al. 2000; Bruscoli, Ferrara & Ciardi 2000; Gnedin and Jaffe 
2000) for direct calculations of the power spectrum of CMB 
variations, C;. It is interesting that results based on simula- 
tions differ in some respects from theoretical expectations. 
Thus, for example, Peebles & Juszkiewicz (1998) demon- 
strate the small impact of moderate nonlinearity which can 
be expected during the period of reionization. The toy model 
discussed in Gruzinov & Hu (1998) shown that the ioniza- 
tion in patches also cannot essentially distort expectations 
obtained under assumption of homogeneous reionization. In 
simulations, however, the noticeable impact of both effects 
was found. The more detailed discussions and comparisons 
of these results with theoretical expectations are needed to 
clarify the validity of theoretical assumptions and represen- 
tativity of simulations. 

Here we turn back to the simplest case of linear theory 
and consider the small scale angular variations of the CMB 
temperature produced by velocity perturbations after the 
reionization of the universe (linear Doppler and Ostriker- 
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Vishniac effects) . Recent discussions of this problem in pub- 
lications cited above use the popular spectral description of 
expected small scale variations of ST/T. This spectrum can- 
not be, however, found with both MAP and Plank missions 
and to do this special extensive observations are required. 
Here we demonstrate that the more simple direct analysis of 
the map of CMB temperature variations provides also valu- 
able information about the reionization epoch and in some 
respect could be more perspective. 

As was shown in Kaiser (1984) and Ostriker & Vish- 
niac (1986), ST/T generated after epoch z = z T where the 
universe last is optically thick, T e (z T ) ~1, are small because 
the different phases of any single plane wave cancel each 
other out, and the main effect is generated at the epoch of 
reionization at z — z ri and/or at redshifts when T e (z T ) ~1. 
For earlier reionization when z ri > z T and the universe is 
optically thick, r e (z ri ) ^>1, the 8T/T averages the pecu- 
liar motion across the Hubble length c/H(z T ) at the epoch 
of last scattering (Sunyaev 1978). In the case, the damp- 
ing of perturbations is defined by the variation of optical 
depth, T e (z), with redshift. For later reionization, z ri < z T , 
the dynamic of reionization becomes important and random 
angular variations of z r i can also be seen. 

The dumping depends upon the velocity coherent 
length, l v , defined by the initial power spectrum. Some char- 
acteristics of density and velocity fields in the models with 
Harrison-Zel'dovich initial power spectrum and Bardeen et 
al. (1986) transfer function were found in Demiahski & 
Doroshkevich (1999; 2001). Here we show that for such mod- 
els this damping is moderate and an amplitude 8T/T ~ 
10" 5 - 3 • 10~ 6 can be expected in minute and subminute 
ranges. 



2 PERIOD OF REIONIZATION 

The processes of reheating and reionization of intcrgalac- 
tic gas have been widely discussed during last years (see, 
e.g., Gnedin & Ostriker 1997; Baltz, Gnedin & Silk 1997; 
Haiman, Rees & Loeb 1997; Tegmark et al. 1997; Abel et 
al. 1998; Silk & Rees 1998; Haehnelt, Natarajan, & Rees 
1998; Miralda-Escude et al. 2000; Haiman, Abel & Rees 
1999; Haiman & Loeb 1999; Shapiro, Iliev & Raga 1999; 
more references in Haiman & Knox 1999) but up to now it 
is far from clarity. 

It is agreed that the reionization is produced by the first 
population of relatively low massive galaxies with 
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(Haiman, Rees & Loeb 1997) where z V i r is redshift of for- 
mation of virialized cloud. The formation of more massive 
galaxies, quasars and black holes is also discussed as at red- 
shifts z ~3 - 5 both bright galaxies and quasars are observed 
(see, e.g., Steidel et al. 1998; Fan et al. 2000). 

The hard problem is a reasonable estimate of luminosity 
of objects responsible for reionization and of the efficiency 
factor of reionization, E, that is a ratio of ionized and col- 
lapsed volumes. Available now estimates of this factor vary 
from E = 114 and up to E = 10 2 [20/(1 + z)f (Haiman & 
Loeb 1999; more references in Knox, Scoccimarro & Dodel- 
son 1998, and Haiman & Knox 1999). Evidently, this factor 



directly relates to the fraction of matter, f m (z ri ), accumu- 
lated by objects with M > Mmin at the redshift of reioniza- 
tion, z ~ z r i, and the mass conservation shows that 



fm{Zri) ~ 1/E(z ri ) 



(2.2) 



Our estimates based on the statistics of high density 
clouds formation (Demiahski & Doroshkevich 2001) confirm 
that for both values of E cited above and for larger mass 
of the heaviest DM particles Mdm >1 GeV the reioniza- 
tion is quite possible at redshifts z r i ~40 - 20 when objects 
with M > M min incorporate ~2 - 8% of the matter. At 
the same redshifts amorphous irregular clouds with a mod- 
erate overdensity above the mean density accumulate ~15 - 
40% of matter. Some part of these clouds can be unstable. 
The latest fraction increases up to ~80% at redshifts z ~4 - 
5 that is consistent with results obtained with high resolu- 
tion simulations (see, e.g., Zhang et al. 1998). The expected 
fraction of homogeneously distributed gas at the period of 
reionization is 



fg{z ri ) ~ 0.4-0.7, 



(2.3) 



for Zri ~20 - 40, respectively. 

These numerical estimates are obtained for the CDM- 
like primordial power spectrum discussed in Sec. 3. They are 
very sensitive to the shape of small scale primordial power 
spectrum and mass of DM particles. More conservative ap- 
proach to the problem of reionization can be found in Ben- 
son et al. (2001) where z r i ~10 is discussed. Observations of 
small scale 8T/T will significantly clarify this problem. 

The reionization and reheating take place at the period 
when only a small fraction of baryonic matter is accumu- 
lated by high density clouds. As was shown in Peebles & 
Juszkiewicz (1998), these clouds do not distort the expected 
variations of CMB temperature. Non the less, the reheating 
of baryonic component suppresses the small scale baryonic 
perturbation in scales lesser then the Jeans scale, Rj, 



Rj 



o s (l 



0.12Vn 

m / - r i 



ft _1 Mpc, (2.4) 



where a s , T4 — T/W 4 K ~ 1 and Q. m are the sound speed, 
the expected temperature of reionized gas and the mean 
density of the universe. It essentially restricts the expected 
variations of the CMB temperature in small scale and, in 
fact, leads to the cutoff in the power spectrum of the CMB 
variations at I ~ 3 ■ 10 4 - 10 5 (Sec. 6). 



3 BASIC COSMOLOGICAL MODEL AND 
CHARACTERISTICS OF THE DENSITY 
AND VELOCITY FIELDS 

3.1 Parameters of cosmological model 

As a basic model, we use the spatially flat ACDM model 
(Jenkins et al. 1998) with parameters 



fir 



0.3, fl b =0.03, ft = 0.7, a v (0) = 650/sm/s, (3.1) 



where £l m & flj, are dimensionless density of dark mat- 
ter (DM) and baryonic components, Qa = 1 — flm, h — 
Ho/WOkm/s/Mpc is dimensionless Hubble constant, and 
a v (0) is a velocity dispersion at redshift 2 = 0. The analy- 
sis of this simulation (Jenkins et al. 1998; Demiahski et al. 
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2001) shown that it reproduces well characteristics ol large 
scale matter distribution observed at small redshifts and, so, 
can be used as a basic model for further investigations. 

Here we are interested in larger redshifts, z >10, where 
the influence of Qa becomes negligible and the main rela- 
tions are simplified. For such redshifts we have 



(3.2) 



a v {z) _ 1.6 ■ 10' 
V3c Vz 




fi m 650km/s 



w 



The velocities of baryonic component suppressed at z ~1000 
are practically restored already at z ~50. 

Distances along a line of sight are measured by the con- 
formal time and linked to the redshifts as follows: 



^^MfA-il + z)- 1 / 2 }, 



(3.3) 
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3.2 Characteristics of density and velocity fields 

To estimate the efficiency of the formation ST/T we use 
two characteristics of peculiar velocity field, namely, the co- 
herence length, l v , and the correlation functions of veloci- 
ties. For the Harrison-Zel'dovich initial power spectrum of 
perturbations, p(k) oc k, and for the Bardeen et al. (1986) 
transfer function, T 2 (x), the velocity coherent length was 
introduced in Demiahski & Doroshkevich (1999) as follows: 



r 2 = 



r = Q m h 



lv ~ p 



kT 2 (k/ko)dk = m-ik 2 ., 



Th 
IMpc' 



1.68p 

Prel 



m_2 



poo 

= / xT 2 (x)dx^ 0.023, 
Jo 



0.023 

m_2 



r'Mpc = 33/i _1 Mpc 



(3.4) 



where p 1 k, p re i are the density of CMB and relativistic 
particles (photons, neutrinos etc.). 

3.2.1 Doppler effect 

The normalized velocity correlation function can be approx- 
imated by an expression (Demiahski & Doroshkevich 1999): 

< Vi(r 1 )v J (r 2 ) > 



F 1 (x)S ij + x i x j F 2 (x), (3.5) 



ri - r 2 . . 

x= — , x = \x\ 

lv 



2 1 r 



dk p(k), 



ji(kl v x) 



kl v x 1 + x + ax 2 ' 



F 2 (x) 



dFi 3_ r 

xdx 2n 2 a 2 J 



where ri and r 2 are unperturbed coordinates of points at 
z — and the parameter a w 0.3 can be expressed through 
other moments of the power spectrum. 

Fits of the perturbation amplitude to the CMB 
anisotropy (Bunn & White 1997) link model parameters 
(3.1) with COBE data and we get: 



„.» 1 3.4„-Mpc(A)(5 f )' 
. Ho<j s 



Q 



15 pK 



wo « 0.65- 



:V3 



where Tq is the amplitude of quadrupole component of 
anisotropy. For the projection of the correlation function 
(3.4) on two directions, 71 and 72, we get: 



F{x,n) = tf v {x) - (1 - ^yllplF^x), 



7172 = P '■ 
Here the function 



I 2 



1-072, !-/*<!. 



(3.6) 



£v(x) = F 1 (x) + x 2 F 2 (x) 



1 — ax 



(3.7) 




(1 +x + ax 2 ) 2 ' 
describe the velocity correlation along a line of sight (for 
/* = !)■ 

To estimates of the variations of CMB temperature at 
the period of reionization we will use the point separation 
measured along a line of sight (3.3) and for the characteristic 
correlation angle we have 



(3.8) 



3.2.2 Ostriker-Vishniac effect 

The cross correlation of velocity and density fields is de- 
scribed by the second order power spectrum of perturba- 
tion as this correlation vanishes for the first order perturba- 
tions. This correlation leads to second order correction for 
the Doppler perturbations and to the Ostriker - Vishniac 
effect. The second order spectrum of baryonic component is 
found in Appendix A. 

The normalized correlation functions, Gi and G 2 , and 
the amplitude of CMB temperature variations for the OV- 
effect, Gb, can be written as follows: 



G ij {y) = (v i {v l )5{v l )v j {v 2 )5{v 2 ))/al-- 
[2G l (y) + y 2 G 2 (y)}8 t3 -y t y,G 2 (y), 



(3.9) 



ri - r 2 2 h Q m 20 



rj « 0.75i?jln ( — ) « Q.YI\ Ta^-— /i _1 Mpc, 

\ajj V \l m Zri 



Gi(y) 



1.2 



(1 + 1,2)0.6' G2 (y)~ (l+J/2)l-6 



2 /0 x 1-2 



= fov«l(0), fov - In (-) (-) (^) 
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Here as before n and r-2 are unperturbed coordinates of 
points at z — and Ob is an amplitude of perturbations. 
The Jeans scale, Rj, were introduced by (2.4). 

For the projection of the correlation function (3.9) on 
two directions, 71 and 72, we get: 

jijjGij = 2 M Gi(z) + (1 -V)^£^G 2 (x), (3.10) 



7172 = /i ~ 1 — # 2 /2, l-/i<l, 
and for the characteristic correlation angle we have 



, ^"jT^-f. (3.11) 

VOJA V 2 " /A 



4 VARIATIONS OF OPTICAL DEPTH 

Now it is generally agreed that the reionization starts with 
formation of separate bubbles number and size of which 
increase with time. At redshifts z = z ri these bubbles 
merge together and the high ionization degree, fi(z r i) ~1, 
is achieved. The maximal optical depth of the universe 
achieved due to reionization can be written as follows: 

n e (x)dx 3/2 



xH(x) 



T Z 



K 



to = i.4- io~ 3 A / — 



3 / 0.3 Mh 



(4.1) 



3/2 j t 
Mx)— > 1. 



Here c & ct c = 6.6 • 10 25 cm 2 are the speed of light and the 
Thompson cross section, the density of electrons n e (z) ~ 
1.2-10" 5 n 6 (l + 2) 3 cm , and a factor Kf takes proper ac- 
count of dynamic of reionization. Thus, Kf — 1 for the in- 
stantaneous reionization at z = z r i and Kf >1 when partial 
ionization is reached at larger redshifts and z 3//2 fi(z) slow 
decreases for z > z ri . In particular, this function integrates 
separate ionized bubbles formed before the full reionization 
and, so, is a random function of line-of-sight. 

Three models of ST/T generation should be considered. 
For earlier reionization when r max >1 we have 



Z T 



zo 



-2/3 



80 



03 



x / 3 /0.02\ 2 / 3 



(4.2) 



In the case, the characteristics of ST/T depend only upon 
parameters of the basic cosmological model. For later reion- 
ization, when Tmax >1 but z ri < zq we have 



zo > z T 



tv-2/3 ^ 



(4.3) 



For r max <1, the generation of ST/T is driven by the dy- 
namic of the process of reionization and final estimates of 
ST/T become more uncertain. In particular, in the case ran- 
dom angular variations of T max can be observed. 

When the CMB radiation is scattered in a cloud with 
a proper size ~ l v and a correlated velocity along a line-of- 
sight v(z), the temperature perturbations of CMB observed 
with angular resolution smaller than a cloud size, 8 c i ~ 10', 
are 

(t 



-J = a c n e (z)l v v(z)/c ; 



(4.4) 



2.3- 10" 



ly 



l+2\ 3 / 2 n b K 2 ^(0) 



33ft- 1 Mpc 



0.015 650/sm/s' 



This amplitude is comparable with estimates of Gruzinov & 
Hu (1998) and with the expected small scale variations of 
the CMB temperature (Sees. 5.1, 5.2). The measurable am- 
plitude of ST/T ~ 10~ 5 can be achieved only within ionized 
bubbles with an extremal velocity <j v ~ 3a„(0). Formation 
of such bubbles is not forbidden even at larger redshifts but 
it becomes more probable in a course of merging of many 
low massive bubbles just before the reionization. For later 
reionization, T max <1, such bubbles could be observed as a 
random angular variation of z r i and Tmax- 

Factors fi(z) and Kf characterize the redshift variation 
of number density of the homogeneously distributed ionized 
gas and must be corrected for probable matter concentration 
within low massive clouds. Small fraction of denser clouds 
provides the reionization but essential fraction of gas can be 
accumulated within irregular clouds with various sizes and 
moderate overdensity above the mean density, 8p/p ~5 - 10. 
As was shown by Peebles & Juszkiewicz (1998), action of 
this effect cannot substantially change the expected optical 
depth, r, and parameters zo and z T but for generality wc 
will include it in our consideration. 

Now the possible formation of such clouds with variety 
of morphology is widely discussed. For example, in simula- 
tions at redshift 2 ~4 (see, e.g., Zhang et al. 1998) up to 
90 - 95% of gas are found to be accumulated within such 
amorphous clouds that explains the suppression of observed 
Gunn-Peterson effect. Efficiency of such concentration at 
larger redshifts depends upon the amplitude and coherent 
length of density field and, so, is sensitive to the mass of DM 
particles and the shape of initial power spectrum. 

Here we will characterize action of this effect by the 
mean fraction of homogeneously distributed gas at redshift 
of reionization, f g — (f g (z)) <1. The expected values of f g 
can vary in wide range, 0.05< f g <1 with redshift, coherent 
length l p and dynamic of reionization. This factor can be 
incorporated in relations (3.1) - (3.4) by the substitution 
of f g Qb in place of Qb that increases our previous estimates 
of zq and decreases estimates of T max . The redshift varia- 
tions of f g (z) can be incorporated in the factor Kf by the 
substitution of f g {z)fi(z)/{f g (z)} in place of fi(z). 



5 VARIATIONS OF CMB TEMPERATURE 

The perturbations of CMB temperature caused by Doppler 
effect are defined as follows: 



(ST 
V T 



>D Jo C 



-Te(z) 



dr e (z). 



(5.1) 



(see, e.g., Hu, Scott & Silk 1994; Peebles & Juszkiewicz 
1998). Here v(z) is a line-of-sight peculiar velocity and, as 
before, r e (z) is the optical depth due to Thompson scatter- 
ing of CMB. For a cloud situated at redshift 2 with a proper 
size 



D cl = cH- 1 (z)Az, Az = H(z)D cl /c « 1, 



(5.2) 



this expression becomes identical to (4.4). In general case it 
defines a random function 8T/T. For Gaussian initial per- 
turbations the probability distribution function of ST/T is 
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also Gaussian and its properties are characterized by the 



For more realistic case of later reionization we have 



variance, a D : 

a D = / 00 dre(z 1 )dr e (^) ^M^ e-^)-^) <r,:i, 
Jo c 



J I U I \ < J v(zi)a v (z2) t . . -T e (zi)-T e (z 2 ) 

3c^ 



Here the function £,,(21,22) describes the correlation of pe- 
culiar velocities along a line-of-sight and, therefore, the ex- 
pected damping of perturbations at the period of reioniza- 
tion and at smaller redshifts. This expression can also be 
obtained directly from Eqs. (25)-(32) of Hu, Scott & Silk 
(1994). 



0.03 



5/6 



(5.10) 



where Kf is given by (4.1), and f g <1 is the mean fraction 
of homogeneously distributed gas at 2 = z ri . The random 
factor Kf characterizes the dynamic of reionization along a 
line-of-sight. 

The measurable variations of ST/T ~ 10~ 5 can be ex- 
pected with a reasonable probability for f g Q. b z r iKf >0.5 
when at least higher peaks in Gaussian field can be re- 
ally observed. The largest clouds ionized at higher redshifts 
magnify Kf and can be seen as random higher amplitude 
patches. 



5.1 Amplitude of CMB fluctuations generated by 
the Doppler effect 

For the basic model, the correlation of velocities along a 
line-of-sight are described by (3.5) with fj, — 1 and 

1 

x — — 




Vo 



1.35 • 10" 



2.0307 

20 n m h y ' 



For kd <SC 1 the expression (5.3) can be rewritten more 
explicitly as follows: 



4^ 



3/2 3/2 

at, « -Tq ~wq I dxidx2^(z T xi,z T X2)ex.p[—To{x 1 ' +x 2 ' ^ 



2.25t^ 3 Wqk d K d , 

> 1, 



(5.5a) 



2 



9 2 2 2 
-T W Z ri 



I' 

Jo 



dXldX2^(ZriXl, Z ri X2)e 



5.2 Amplitude of CMB fluctuations generated by 
the Ostriker-Vishniac effect 

Amplitude of the CMB fluctuations generated by the 
Ostriker-Vishniac effect can be found in the same manner: 



2 _ 2 KjKjf ov 

<Jov — o D — — — , 

kdKdz^ 



(5.11) 



KjKj 



dz\dz2G\(z\, 22) -t( Z1 )-t(z 2 ) 
(l + 2l)(l + 2 2 ) 6 



KJ = 



<JOV 



rj 



no 



0.7- 10" 



T4 2tT' 



Kj « 53. 



n9 ( h n m \ 15 /2o\ 1/4 »o_ 



10- 3 

Here the function G\ (21,22) describes the correla- 
tion of perturbations along a line-of-sight, parameters 
Vo,r.j, fov, kd,Ko&o?d were introduced in (3.3), (3.9), 
(5.4), (5.5) and (5.10). 



2 . 25tq Wq z^.i kd Kd , 



< 1. 



(5.56) 



Here kd describes the main damping of perturbations while 
the factor Kd <1 takes proper account of damping and 
weights it along a line-of-sight. Kd weakly depends upon 
small redshifts, slowly decreases with kd and Kd «0.25 - 

0. 5 for the range of interest, kd ~0.01 - 0.1. For kd <0.01 
Kd decreases oc kd rn(l/Ko). Parameters wo, and r were 
introduced by (3.2) and (4.1). 

For the hypothetical case of earlier reionization, T max ^> 

1, 2 r = 20, we have 



n ~3 m- 5 f f ^ b Y /6 ( °- 15 ^ 

aD ~ 3 ' 10 I O03 J lo^J 



o-.(O) 
650km/ s ' 



(5.9) 



This result confirms that for so earlier reionization, z ri > 
z T ~ 80(f g n b h/0.02y 2/3 (Qm/0.3) 1/3 , r max >1 the mea- 
surable variations of ST/T ~ 10~ 5 can be expected with a 
high probability even if the velocities of baryonic component 
are not yet restored and cr„(0) <650km/s. In the case, the 
amplitude depends mainly upon basic parameters of cosmo- 
logical model and random angular variations of z ri caused 
by the random spatial distribution of earlier ionized bubbles 
are not seen. 



6 ANGULAR CORRELATION FUNCTION 
AND POWER SPECTRUM OF CMB 
VARIATIONS 

6.1 Power spectra 

The power spectra of CMB variations can be found directly 
from the power spectra of perturbations (Vishniac 1987, Ap- 
pendix B & C): 



^^6.4.10- 12 ^T 2 (1V 
2?r a \aj ' 



I 



2tt 



(6.1) 



(6.2) 



for Doppler and Ostriker-Vishniac effects, respectively. In 
both cases the power spectra of CMB variations reproduce 
well the transfer functions of perturbations, namely T 2 &cT b . 
The transfer function of baryonic perturbations, T b , is intro- 
duced in (A. 8). As was shown in Appendix B & C for larger 
I, Ci oc z r i for the Doppler and Ci oc yfzr~i for the Ostriker- 
Vishniac effects, that verifies the domination of period of 
reionization for both effects (Vishniac 1987). 

These power spectra are plotted in Fig. 1 together with 
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Figure 1. The spectrum of the CMB temperature anisotropy 
generated by the Dopplcr (dots) and OV-effects for ctj = 
0.05 (triangles) & 0.01 (squares) as compared with the expected 
primary CMB anisotropy (solid line). 



the spectrum of the expected primary CMB anisotropy. The 
variations of the amplitude with cosmological parameters 
and redshifts of reionization are described by (5.10), (5.11), 
(B.4) and (C.3). 



6.2 Correlation functions 

The angular correlation function of the CMB fluctua- 
tions generated by both effects, w(0), can be found from 
the general expressions (3.6) and (3.9) using functions 
£u,-p2, Gi & G2 and the point separation 

2 




/a 1 



/a 1 



Zl 



Z2 



(6.3) 



for the Doppler effect and similar relation for Ostriker - Vish- 
niac effect. 

The normalized correlation functions are plotted in 
Figs. 2 & 3 for z ri = 20 and 40. For 40 > z ri > 10 these 
functions can be well fitted by expressions 



W(X) « (1 + x 2 ) 1 
7 d « 0.43, 6>o « 



4.20„ 



9/0o, 
« 30', 



(6.4) 



7 j « 0.15, O « 40J « 8", 
for Doppler and Ostriker- Vishniac effects, respectively. 

6.3 Characteristics of peaks 

The small expected amplitude of the CMB variations (5.10) 
and (5.11) indicates that the observational determination of 
the power spectrum Ci is hard problem and at the first step 
only high peaks with ST /T ~ 2 - 3 a Q v ~ 10" s could be 



Figure 2. Correlation functions of CMB fluctuations generated 
by the Doppler effect for z r i = 20 (points) and z r i = 40 (stars) 
together with fits (6.4). 
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Figure 3. Correlation functions of CMB fluctuations generated 
by the Ostriker- Vishniac effect for z r i = 20 together with fits 
(6.4). 



actually observed. This means that the peak distribution is 
also an important characteristic which helps to discriminate 
the cosmological patches from the noise. 

The expected shape of peaks is described by the correla- 
tion function (6.4) whereas their number density is given by 
the Minkowski Functional. Application of this approach to 
the variations of CMB temperature have been discussed in 
Winitzki & Kosowsky (1997), Schmalzing & Gorski (1997), 
Novikov, Feldman & Shandarin (1998), Dolgov et al. (1999). 
For the Gaussian distribution of the temperature variations 
with the correlation function (6.4) we get for the number 
density of positive peaks higher then threshold amplitude, 
is = ST/T/ar, 



(n ma x(> v)) = 



1 + 7 
2tt^0 ( 2 



/("), 



(6.5) 
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where the function f(u) characterizes a threshold depen- 
dence of Umax- This dependence helps to select and iden- 
tify the peaks under consideration for larger statistics of the 
peaks with different v (Novikov et al. 2001). 

The expected number density of higher Doppler peaks 



{n max (> v = 1)) « (113')" 
{n max (> v = 2)) « (250')" 



(6.6a) 
(6.66) 



The expected distribution function of CMB fluctuations 
generated by the nonlinear OV - effect is also close (but not 
identical) to the Gaussian one. Indeed, as is well known, the 
distribution function for the product of two normal variates 
is 

/(i/) = 1/tt Jfo(M) oc arp(-M), \v\ » 1 

where Ko is the modified Bessel function (Stuart & Ord 
1994). However, the integration of perturbations along a 
line-of-sight makes again the distribution of CMB fluctu- 
ations similar to the Gaussian one. The expected deviations 
from the Gaussianity are small for \u\ <1 but hold more 
strong for larger The degree of these deviations depends 
upon the actual number of "clouds" within the integral (5.1) 
along a line-of-sight and, therefore, upon the redshift of 
reionization. It increases for later reionization when small 
amplitude of the CMB variations is also expected. 

The application of Minkowski functional to some non- 
Gaussian statistics of CMB had been illustrated in many 
publications (see, e.g., Coles & Barrow 1987, Coles 1988). 
Their results verify that, for example, deviations between 
the Gaussian and \n statistics with n > 30 do not exceed 
10 - 12%. The significant difference between characteristic 
angles as given by (3.11) and (6.4) (6o ~ 49j) as well as a 
large difference between aov (5.11) and the possible con- 
tribution of a separate cloud with a size ~ r.j indicate that 
the number of "clouds" along a line of sight is 3> 20 and the 
expected degree of the deviations from the Gaussianity can 
be moderate even for the highest peaks. 

Because of this we will apply expressions (6.5) for the 
first estimates of surface density of higher OV-peaks as well. 
We have: 



{n max (> v = 1)) m (37")" 
{n max (> v = 2)) « (82")" 
{n max (> v = 3)) « (5')~ 



(6.7a) 
(6.76) 



(6.7c) 



The same expressions describe the mean number density of 
the negative peaks, n min (< v = —1), n min (< v = —2) 
and n m in{< v = —3). These estimates show that with 
the angular resolution ~ 8 — 10" peaks with the amplitude 
8T/T ~ 10~ 5 can be observed with a reasonable probability. 



6.4 Negative intensity patches 

The discrimination of the cosmological 8T/T becomes more 
easy if we concentrate more attention on the investigation of 
negative intensity patches (Dubrovich 2000, Springel, White 
& Hernquist 2000). The expected Gaussian distribution of 
cosmological 8T/T predicts equally probable patches with 



both positive and negative intensity, but the negative in- 
tensity patches cannot be produced by any background or 
foreground sources excluding only Sunyaev-Zel'dovich (SZ) 
effect in clusters of galaxies some of which can be identified 
with their x-ray radiation (at least at z <1, Refregier et al. 
2000). 

To discriminate the noise, SZ and Doppler patches their 
different frequency dependence can be used (Sunyaev & 
Zel'dovich 1980). The Doppler shift discussed above does 
not change the spectrum of CMB and the frequency varia- 
tions of intensity of Doppler and OV patches are 



SF D ( X ) = F 



4 _ Y 

X e x 



5T 



(1 - e-x)2 T ' 



X = hv/kTj, 



(6.8) 



Fq = 2.7 ■ 10 15 ergs(cm 2 s sr 



Hz) 



T-, = 2.73F 



where \ ls dimensionless frequency of radiation. The max- 
imal intensity of these patches, 8F/F « A.9(ST/T), is 
achieved at \ ~3.8. It drops up to SF/F ^ ST/T at \ ~1 
and x ~8.6 . 

The frequency dependence of SZ effect (Zel'dovich & 
Sunyaev 1969; Weller, Battye & Albrecht 1999) is 



SF az (x)=F y c 



4 - 

X e 



(l-e-x)2 



1 + e~ x 
1 - e-x 



(6.9) 



where y c is the dimensionless Compton parameter. This in- 
tensity is negative at \ <3.8 and positive at x >3.8, max- 
imal SF sz /Fo = 6.8y c and minimal 5F SZ /F — —4.1y c are 
achieved at x =6.5 and \ =2.3, respectively. 

If observations are performed with a set of wavelengths 
in range 1 cm> A >0.1 cm then the different frequency 
dependence of (6.8) and (6.9) allows to discriminate between 
patches produced by Sunyaev-Zel'dovich and Doppler and 
OV effects. 

The free- free emission of reionized plasma can, in prin- 
ciple, restrict the existence of negative intensity patches. It 
depends upon many factors such as the temperature, ion- 
ization degree and small scale clustering of plasma at high 
redshifts estimates of which are strongly model dependent. 
However, the available estimates of the observed intensity of 
the free-free emission 

Fff < 3 • 10~ 19 ergs(cm s sr Hz)' 1 = 10~ 4 F (6.8) 

(Bartlett & Stebbins 1991) confirm that within discussed 
range of wavelengths the cosmological patches can be ob- 
served. More detailed discussion of the background and fore- 
ground can be found in White et al. (1999), Cooray, Hu & 
Tegmark (2000), and Refregier, Spergel & Herbig (2000). 



7 SUMMARY AND DISCUSSION 

Different methods of investigation of the "dark age" in the 
history of the universe evolution between z =1000 and z = 5 
were discussed during last years (see references in Sec. 1) 
but most attention has been concentrated on the traditional 
spectral approach to the investigation of small scale varia- 
tions of CMB temperature that requires extensive observa- 
tions with high precision. The expected amplitudes of cos- 
mological ST/T obtained in Sec. 6 show that the more simple 
direct observation of higher positive and especially negative 
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intensity patches can successfully complement the spectral 
investigations. 

The patches can be observed with already available 
techniques using the aperture synthesis and/or interferome- 
try (White et al. 1999; Dubrovich & Partridge 2000). Radio- 
frequency interferometry offers the advantage of virtually 
complete freedom from atmospheric emission and other sys- 
tematic effects. Now such telescopes as the ACTA (Aus- 
tralia) and VLA (USA) in cm as well as IRAM (Plateau de 
Bure) and JCMT (Hawaii) by SCUBA in mm and submm 
ranges can be used to search for the patches. 

To discriminate reliably between Doppler, OV and SZ 
patches, wavelength A <2 mm should be used. In future 
most promising will be ALMA, to be completed in the next 
decade. It will operate in all the atmospheric windows from 
1 cm to 0.45 mm and will be able to detect sources at A <0.1 
cm with the angular resolution and sensitivity of the present 
VLA at cm wavelengths. 

Higher amplitudes (5.9), <jt ~ 3- 10 , expected for the 
hypothetical case of very early reionization, z r i > Zo ~80, 
depend mainly upon cosmological parameters, Q. m h 2 and 
cr„(0). Perhaps, so early reionization is possible for the tilted 
power spectra with an excess of power in small scales. In 
the case, distortions of the first generation of ST /T can be 
essential even at angular scale ~ 1° (Griffiths, Barbosa & 
Liddle 1999; Tegmark & Zaldarriaga 2000). 

More realistic scenario of later reionization, z r i ~40 - 
20, predicts the existence of patches with or ~ 0.5 • 10~ 5 as 
given by (5.10), (5.11) which can be observed at least as a 
system of the separate highest peaks. The amplitude is larger 
for earlier period of reionization due to growth of both z r i 
and fg(z r i). Reionization at z r i >20 is possible for the tilted 
power spectra with an excess of power in small scales and/or 
for the Harrison- Zel'dovich power spectrum with extremely 
heavy DM particles, Mum ~ 10 18 - 10 21 keV. Existence of 
such particles is now discussed as a possible interpretation 
of observed high energy cosmic rays (Bertou, Boratav, & 
Letessier-Selvon 2000; Nagano & Watson 2000). In the case, 
possible random angular variations of z r i caused by random 
merging of earlier ionized bubbles can also be observed as 
the separate highest peaks. 

The angular distribution of the patches is driven by 
the spatial correlation of the peculiar velocity field. The 
same field is also responsible for a large scale distribution of 
galaxies and their strong concentration within filaments and 
walls (see, e.g. Demianski & Doroshkevich 1999, Demianski 
et al. 2000). This concentration is observed in deep surveys 
such as the Las Campanas Redshift Survey (Shectman et 
al. 1996) and Durham/UKST redshift Survey (Ratcliffe et 
al. 1996) and is well reproduced in numerous simulations. 
For example, the mean separation of richer observed walls, 
D sep ~ 50 — 60/i _1 Mpc, is about twice of the velocity co- 
herent length, l v , as defined by (3.4). This implies that the 
correlated distribution of negative patches and, perhaps, the 
appearance of the network of patches with few angular min- 
utes sizes can also be expected. The correlation of higher 
peaks was discussed in Heavens & Sheth (1999). 

Similar network of negative patches with the angular 
sizes 9 ~1 — 2' and angular separation 8 3ep ~ 5' was firstly 
observed by Carlstrom, Joy & Grego (1996) with 30 GHz re- 
ceiver of Owens Valley Radio Observatory in a course of in- 
vestigations of Sunyaev-Zel'dovich effect toward two clusters 



of galaxies. The statistical reliability of these small ampli- 
tude negative patches is not enough, and these observations 
should be repeated with better sensitivity and in wider range 
of wavelengths, 1 cm> A >0.1 cm. The frequency depen- 
dence of the amplitude and statistical properties discussed 
in Sec. 6 allows to test the attribution of the patches as 
cosmological and generated at redshifts z >10. 
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where Vi and 8 are peculiar velocities and relative density 
perturbations. Here 

$y(fc) = / d 3 q{Vi(k)S(k)v 3 {q)8(q)) = (Al) 



a«-^)*i(*) + ^*2(*). 



*a(*) 



_p(g)p(|k-q|)-^ 



(k- 



(kq) (kq)\ 



The function $1 describes the Vishniac - Ostriker effect 
while the function $2 gives the second order contribution 
to the Doppler effect. Both functions characterize the DM 
perturbations and at k — > $1 oc $2 oc k 2 while for k — > 
00, $1 oc $2 oc k~ z and they are sensitive to the small 
scale nonlinear damping of perturbations. These functions 
are well fitted by expressions: 



oc x 2 (l + 2.5x + Ax 2 + 0.37a; 2 - 5 ; 
$ 2 {x) oc x 2 (l + 1.45a; + 5.37i- 2 + 1.25a; 2 



(A2) 
(A3) 



where x = k/ko. The function fc 2 $i(fc/fco) plays a role of 
the power spectra for OV-effect. 

To describe perturbations of baryonic component and 
to take into account the pressure of reionized gas we will 
substitute in equations (AT) instead of p(k) 

p(k)j' 2 (k Rj ), Rj = ^L£l M ^Ifc-iMpc. (A4) 

Here a s and T4 = T/W K ~ 1 are the sound speed and 
the expected temperature of reionized gas. This factor sup- 
presses the small scale perturbation and instead of (A. 2) we 
get for perturbations of baryonic component 



3>b(x) ~ 5>i(:r)exp(— O.GctjX ), 



(A5) 



x = k/k , aj = k Rj = 10 Tr^A/Ti--— — • 

0.7 V 0.3 Zri 

With this spectrum we get for the normalized spatial 
correlation functions, Gi and G2, the transfer function Tb 
and the amplitude of OV-effect, at: 



Appendix A 



Gij(y) = (vi(ri)<J(ri)vj(r2)5(r 2 ))/CT6 



(A6) 



Second order power spectrum [2Gi(y) + y 2 G2{y)\5ij - yiVjG^y), y = — = 



for the Ostriker- Vishniac effect 

The cross correlation of velocity and density fields is 
described by the second order power spectrum of perturba- 
tion as this correlation vanishes for the first order perturba- 
tions. This correlation leads to second order correction for 

the Doppler perturbations and to the Vishniac - Ostriker 3 r°° ji(kri2) 

effect. Gl (y) = — dxx 2 T b {x) 



The second order spectrum can be written as follows: 



rj rj 



r., w 0.75J?jln (— ) w 0.17, T A ^- — /i _1 Mpc 

\OLjJ V "m Zri 



m b J v ' kn 2 (1 + t/ 2 ) - 6 



(^(nJ^nJw^raWra)) = / d 3 fc$ tJ (fc) exp[ik(n - r 2 )], G 2 {y) = -—( dxx A T b (x) j2 ' '' ' ,j1 - 



m b J Q ° v ' (fcr 12 ) 2 " (1 + y 2 ) 1 - 6 
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To Or) 



x 2 exp(— 0.6a 2 jX 2 ) 
(1 + 2.5a; + 4a; 2 + 0.37a; 2 - 5 ) 2 ' 



x = k/k , (A 7) 



m b 



2 



= / dxx 2 T b « 0.1161n 2 (a. 7 ), 
Jo 

= / d 3 k$ b (k) = fov<£, 
Jo 



(A.8) 
(A.9) 



to 



3.7- W^aj^/m; 



, ln /J_\ (h_\ 2 f^ m \ 12 T Q 
' n \aj) 1 0.7/ V0.3/ 15/j.K' 



Here as before ri and r 2 are unperturbed coordinates of 
points at z = and a b is an amplitude of perturbations. 

Appendix B 

Power spectrum of CMB 
generated by Doppler effect after reionization. 

To obtain the power spectrum of CMB perturbations 
generated by the Doppler effect we use the general expres- 
sion for the angular correlation function 



dz 1 dz 2 (^ v -(l-v 2 )'^F2), (B.l) 



where r)\ = 77(21), 772 = rj(z2) and structure functions £„ and 
F2 are expressed through the power spectrum of primordial 
perturbations: 



rri-Q 



Fo = 



xT 2 (x)dx[j (r,)~2j 1 (r l )/r l ] 



poo 

■ / x 3 T 2 (x)dxj 2 ( V )/ V 2 , (B.2) 
Jo 



V = \f)(zi) - 77(22)!, 77(2) w (M ( /a - 



a; = fc/fco, a = 



ft /n 



2c ^ 2 .io 3 - ,. . 

HoV^Z 0.7 V 0.3 




and T 2 (x) is the transfer function of primordial perturba- 
tions. 

With the 'summation theorem' (Gradshteyn & Ryzhik 
1994, 8.533), we get for the power spectrum of CMB pertur- 
bations: 

r 3(i 



ft 



[ dxxT 2 (x)f?(x), (B.3) 

In, (r,\ 

, J7 = ?7(2) 



dr\ 



This expression is identical to that obtained by Vishniac 
(1987). It can be applied to any primordial power spec- 
trum and any history of reionization and indicates that Ci 
depends upon both these factors. We can consider fi = 
exp[—r(z)]dz / dr) as a slowly varying function and rewrite 
(B.3) as 



1 ^Zyi I 



f°° dx T 2 ( x\ 
Jo x 



6z r 



l 2 Cl _ 3 Zri l_ T 2 fj_\ 

27T 7T ttm-2 a \aj 



o. 5 It 2 (L) 



zh_(17 0.3 



(5.4) 



Appendix C 

Power spectrum of CMB generated by 
the Ostriker- Vishniac effect after reionization. 

To obtain the power spectrum of CMB perturbations 
generated by the Ostriker- Vishniac effect we use the general 
expression for the normalized angular correlation function 



C - 



dz\ dz2 
1 + 21 1 + 22 



(CM) 



where 771 = 77(21), 772 = 77(22) and structure functions G\ and 
G 2 are expressed through the power spectrum of primordial 
perturbations (Appendix A) With the 'summation theorem' 
(Gradshteyn & Ryzhik 1994, 8.533), we get for the power 
spectrum of CMB perturbations: 



C, = 



gi 



6T 
m b 



dxx 2 T b (x)gf(x), 



(C.2) 



dz 



r(z) 



jl (v) 



,o - ■ ~ V 
This expression is identical to that obtained by Vishniac 
(1987). The analytical expressions for G can be found for 
T b (x) cc x°,x~ 2 ,x~ 4 , ... . They show that the approximate 
general expression for G can be written as follows: 



l 2 G 
2tt 



18 Jz, 



7T m b 



(C.3) 



7 (£) n (I) v t? (1^7) 



Numerical factor ~1.5 in (C.3) is found by numerical inte- 
gration of (C.2). 



© 0000 RAS, MNRAS 000, 000-000 



